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Abstract: We investigate the discretization errors affecting correlators of the energy- momentum 
tensor T^jy at finite temperature in SU(A'^c) gauge theory with the Wilson action and two different 
discretizations of T^jy. We do so by using lattice perturbation theory and non-perturbative Monte- 
Carlo simulations. These correlators, which are functions of Euclidean time xq and spatial momentum 
p, are the starting point for a lattice study of the transport properties of the gluon plasma. We find 
that the correlator of the energy j (I^xTqq has much larger discretization errors than the correlator 
of momentum J d^xTok- Secondly, the shear and diagonal stress correlators (T12 and Tkk) require 
Nt > 8 for the Txq = ^ point to be in the scaling region and the cutoff effect to be less than 10%. We 
then show that their discretization errors on an anisotropic lattice with aa/a-r = 2 are comparable to 
those on the isotropic lattice with the same temporal lattice spacing. Finally, we also study finite p 
correlators. 
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1. Introduction 

Energy-momentum tensor (EMT) correlators in finite-temperature QCD provide a way to study the 
spatial correlation of fluctuations in energy density, pressure and entropy [1]. They thus allow us to gain 
insight into the structure of the quark-gluon plasma (QGP), thereby going beyond its thermodynamic 
properties. The time-dependent EMT correlators are related to the transport properties of the QGP. 
Phenomenological upper bounds for the shear viscosity to entropy density ratio, Tj/s < (see [2] 
and Refs. therein), are derived by comparing hydrodynamic calculations of elliptic flow to heavy ion 
collision data. This result suggests the picture of a strongly coupled plasma around 2rc, where Tc is 
the QCD crossover temperature. A significant effort is underway to constrain the transport properties 
of the gluonic sector non-perturbativcly from first principles [3, 4], and given the high computational 
cost of determining the EMT correlators in Monte-Carlo simulations, it is important to optimize the 
choice of lattice action and discretization of the EMT. This is the subject of this paper. 

The gluonic correlators have been computed at treelevel in the continuum [5]. Here we compute 
the correlators at finite lattice spacing using the (anisotropic) Wilson plaquette action [6] and two 



discretizations of the EMT. The strategy is to determine the kinematic regime where the treelevel 
discretization errors are below the 10% level. We then expect treelevel improvement to further reduce 
the discretization errors to the few percent level. This is a reasonable target, given the statistical 
accuracy being currently achieved in lattice simulations [4]. 

The outline of this paper is as follows. In section 2 we give the definitions of the lattice actions 
and discretizations of the EMT, and derive the treelevel formulas for the two-point functions on the 
lattice. Section 3 is devoted to analyzing the treelevel cutoff effects for various correlators. Section 4 
describes how the treelevel discretization errors can be removed from non-perturbative Monte-Carlo 
data, and in section 5 we discuss to what extent treelevel improvement is successful. We end with some 
concluding remarks. 

2. Definitions and master formulcis 
2.1 In the continuum 

The continuum Euclidean energy-momentum tensor for SU(A^c) gauge theories reads 



C/ii/j/iv (iio summation) will be sometimes abbreviated as C^u, and Cqq is the correlator of the trace 
anomaly 9, normalized as in Eq. (2.5). In the continuum, these correlators were calculated to leading 
order in perturbation theory in [5]. The extent of the time direction is denoted by Lq = l/T. 

2.2 On the lattice 

On the anisotropic lattice with spatial lattice spacing and temporal lattice spacing a,-, the Wilson 
action reads 




(2.1) 
(2.2) 
(2.3) 
(2.4) 



We will study the dimensionless, finite-temperature T Euclidean correlators 



c^,,,pa{xQ, p) = r-5 / d^x e'p"' (r^,(xo, x) r^,(o)) . 



(2.5) 




(2.6) 



X 



In numerical practice, it is convenient to parametrize these parameters as 




/3r = /3?0. 



(2.7) 



At treelevel, = o-ajo-T- We use the notation 




SkU Sr = y2^0k, 5^i.(x) = 4:ReTr{l - P^i.(x)}, 



(2.8) 



k<l k 



(2.9) 



The lattice spacing in the four directions are denoted by a^u in order to maintain the symmetry 
among different directions as long as possible. At the end of the calculation we will set oi = 02 = 



03 = Oo- and uq = Ur- The vectors fl are defined as unit vectors along the lattice axes. We employ the 
summation convention for color indices, but not for space-time indices. We use the standard notations 

d^f{x) = -{fix + a^A) - fix)), d;fix) = -ifix) - fix - a^fi)), = + d;) , (2.10) 

We introduce the perturbative fields A^ix) by 

?7^(x) =e^o"''^''(''). (2.12) 

Using an antihermitian set of generators normalized by Tr {T^T^} = — we define A^ix) = A^(x)T'^. 
The latter has the covariant-gaugc propagator 

{A;{x)Atiy))o = S^'J^^/- J, {<^.^-(l-Ao^)^|- (2-13) 

The expectation value (. . .)o is taken in the free theory. Here and in the following, (O1O2) (whether in 
the free theory or not) will always be understood to be the connected two-point function. The Brillouin 
zone is ;B = 'Si/j.B^ with = [— vr/a^, vr/o^]. The lattice field strength is defined by 

P/xuix) = explgoa^a^Fi^^ix)}. (2.14) 

To leading order, one has 

i^M-(x) = di,A,ix) - d^A^ix) + 0(50). (2.15) 

2.3 Clover discretization 

We introduce Q^uix) and F^^ix) as in [7], 

Qiivix) = P^uix) + PfMi^ix - afj^fi) + Pnuix - a„v) + P^^ix - a^^/i - a„v), (2.16) 
%vix) = -^iQ^,uix) - Q.^(x)). (2.17) 

We define F^^ = goF^^T'', implying for instance -4TV {i?,^(a;)F^^(a;)} = F^Ax)F^u{x). Then 

"0 

K'^ix) = \ ^Ml{x) + Alix - v)) - dMlix) + Alix - ft)) + 0(50) 

Using Wick's theorem, one finds^ (d^ = — 1) 

( F^,ix)F^Ax) F'^0iy)F^5(y) )o 

{F^uix)Fapiy))o{Fpaix)F^siy))o + {F,,Ax)F^5iy))o{Fpaix)Fapiy))o 



U(l) 



^If Of are linear combinations of the gauge fields, Of = Xi,aAf„ then by Wick's theorem 

53((0?Of ) (O!oi)>o,conn = 53(0?0?>o{02"02)o + {Of02)o{0^0^)o 

a,b a 
= dA((Ol03){0204) + (Ol04>(0203))u(l)- 

The second equality follows from the fact that (Of 02)0 is independent of a and equal to the corresponding correlation 
function in the U(l) gauge theory. 



Then, in Feynman gauge Aq = 1, we obtain 

{%v{x)Fap{y)) = dA <^,j,ua(}{x - y) , 

(l>Huap{x) = 5^pfl^{x) + S^ocf^pix) - (5^/3 //fa (x) - S^af^pix) , 



JB 



COS {puay/2) Pa Pn- 



(2.18) 
(2.19) 

(2.20) 



'^'^^ ' - Js (27r)4 p2 
Finally, 

{F^u{x)FpAx) F^p{y)Fjsiy))o = dA ^tivap{x - y)(t)pajs{x -y)+ (t>nu'i5{x - y)(t)paai3{x - y) . (2.21) 

This correlator is gauge-invariant and therefore independent of Aq. To go over to mixed propagators 
(which are functions of (xq, p)), we introduce the spatial Fourier transform of 4'ixvap{x), ^'nuafiixo, p) = 
alJ2^<l>^^'^ap{x) e^P-. Then 



alj^e^^y { F-Am^M F^;3(xo,y)Fi(xo,y) )o 

d3p 



(2.22) 



y 

= dA 



(2vr)= 



(t'piuaf){xo,Y*)(t>pcT-i5{xQ, -(P + q)) + 4>^lu^5{xQ,V)4>paaj3{xQ, "(p + q)) 



with = Bi X B2 X B3. Eq. (2.22) is the master formula from which we will derive all results in 
section 2. 

Explicitly, in the clover case, 



, , , , dpo e^fo^o 



Bo 



27r j5g + p2 



(2.23) 



cos^(a;^p,y/2) pp, Pa + Spa cos^{apPp/2) p^ p/s 
Spf3 cos^ {apPp/i) 

Pv Pa "i/a 

cos {a^pvll) Pp pp 



At finite temperature, one is to replace the po-integral in ^ by a Matsubara sum, 

dpo 
IB, 27r 



/ 

J Br 



Pq = 2'KnT, n & 7i. 



(2.24) 



It is worth noting that the scalar propagator in the mixed (xo,p) representation can be calculated 
explicitly, even at finite lattice spacing [8]: 

ar cosh(Lj(iLo - xq)) 



1 



E 



o«poa:o 



Lq ^ po + w^ 2 sinh(a)aT-) sinh(|c2;Lo) 



■ sinh 



(2.25) 



2.4 Plaquette discretization 

The gauge-invariant continuum operator \F^i,{x)F'^v{xy\coxA can be discretized gauge invariantly as 



Re TV {1 - Pp,{x)} = F;,(x)F;,(x) + 0(a^,J , 



(2.26) 



by the lattice field F^^x), where Fp^x) = F^^(a;)r". Then Eq. (2.22) still applies if all the F are 
replaced by F's, provided fpai^) (^^^ ^Q- 2.20) is replaced by 

J4„ 



fpaix) 



r d'^p e'P^ 



Pa Pp- 



(2.27) 



3. Investigating cutoff effects 



In this section we investigate the numerical size of cutoff effects for two different discretization schemes 
of the energy momentum tensor correlators. The length of the temporal direction is set equal to 
Lq = 1/T, where T is the temperature. There is no finite temperature at which perturbation theory 
correctly describes finite-volume effects, since those are related to magnetic screening, which is a non- 
perturbative effect. For that reason, the spatial volume is kept infinite. 

We use Eq. (2.25) and perform the spatial momentum integral in Eq. (2.22) by Gaussian quadrature 
with a target relative accuracy of 10~^. We checked a sample of the results against an extrapolation 
of the finite-volume momentum sums, and by comparing the latter sums to Monte-Carlo data at very 
high j3 value. 

3.1 Correlators of conserved charges 

We start with the isotropic lattice. Figure 1 displays the treelevel lattice correlator of the energy 
operator, ^oo(a;). The trace anomaly is formally 0{as) and does not play a role at this leading 

order. Since in the continuum 

C^oo,oo(xo, 0) = ^ / (Too(x)Too(O)) = ^ Vxq / 0, (3.1) 

the departure of this correlator from a constant Cy/T^ = 47r^dA/15 is a measure of discretization errors. 
The discretization errors fall below the 10% level only for xo/a > 5 for Nt = 20. The reason for this 
large cutoff effect is that generically Cf^iy_pcr{xo,p) fall off as x^^ , and this singularity must cancel in 
the case of f d'^xTooix). It is not surprising that this large cancellation only takes place for xq 3> a, 
where time-translation is effectively restored as a symmetry. 
For the momentum density operator, 

I Yl CokM^o, 0) = ^ E / (rofc(x)rofe(o)) = ^ \/xo^ o, (3.2) 

k k 

the situation is significantly better (see Fig. 1, bottom panel), although the same cancellation has to 
take place. At N^- = 20, discretization errors are well below 10% for xq > 4a. 

3.2 The tensor channel = 1) 

In the continuum and infinite spatial volume limit, the following equality holds by rotational invariance: 

^ J Sx ((Tn - r22)(0) (Tn - T22){x)) = j d'x {Ti2{0)Ti2{x)) , Vxq. (3.3) 

In finite spatial volume, the two correlators differ even in the continuum. On the infinite cubic lattice, 
discretizing either side of Eq. 3.3 yields a correlator that approaches the continuum limit with different 
O(a^) discretization errors. Figure 2 shows a comparison of the discretization errors affecting these 
two schemes. With the clover discretization, at a given value of Nt, discretizing T12 yields smaller 
discretization errors than discretizing ^(Tn — T22), for all values of xq. 

With the plaquette discretization, Eq. 2.26, one is only able to treat the diagonal elements of T^n, 
and in general it leads to significantly larger discretization errors than the clover discretization. If E 
refers to the chromo-electric field and B to the magnetic field, the EE and BB terms in the correlator 
are defined at integer values of x^jaT-, while the EB term is defined at half- integer values of x^jaT- 



For a function falling off as Xq^, without a careful treatment this mismatch leads to large O(a^) cutoff 
effects. An appropriate scheme [9] is to compute separately Cbb, Cee at integer values of xq/q and 
Ceb at half-integer values. One can then obtain the two-point function of ^(Tii —T22) by interpolation, 
which is treated as part of the treelevel improvement scheme (to be discussed in section 4). However, 
at treelevel, Ceb vanishes at the midpoint xq = Lq/2, and Cbb = Cee, so it is straightforward to 
compare the plaquette discretization scheme to the others at that point. As figure 2 shows, its cutoff 
effects are almost identical to those obtained with the clover discretization of T12. 

In all cases, > 8 is necessary for the cutoff effects to be less than 10% and in the O(a^) scaling 
region. 

3.3 The scalar channel (^ = 1) 

We consider two ways to evaluate the zero-momentum two-point function of ^j^. T^jt on the lattice. 
They differ, for large N-j- and xo/ar, by O(a^) terms. 

1. The first way then consists in discretizing ^ Ca^kk directly. Using the notation introduced 
earlier, one may rewrite it identically as 

j d^x{Tuio)Tkk{x)) = j d'xUe^o{^)eo^{x)) - ^{e{{))eoo{x)) + ^{e{{))e{x))\ . (3.4) 

i,k 



2. Alternatively one can exploit the conservation of energy to write 

s — 



Y <^iiM^o, 0) = ^^7^ + Ceeixo, 0) (xo > 0). (3.5) 



i,k 

We have used the standard expressions for entropy and specific heat, s = ^ and = 
Here one computes the trace-anomaly correlator in a given discretization, and subtracts the 
thermodynamic function appearing on the right-hand side, either at the same value of Nj- or 
already extrapolated to the continuum. 

3.3.1 Asymptotic temperatures 

At high temperatures, Cy ~ 3s and we have 

6s — Cy An^dA 



1-— — + 0(a|) 

4 TT 



(3.6) 



r3 15 

In particular, this quantity is positive and straightforward to compute non-perturbatively. Since the 
two-point function is formally 0(0:5), ^1- i^-^) implies that the leading expression for Ca^kk is simply 
Eq. (3.6), independent of xq. The correlators in the scalar channel thus have a large xo-indcpendent 
contribution. At high temperatures, the choice between strategy (1.) and (2.) amounts to deciding 
which of T^(6s — c,,) or f d'^x{9QQ(0) Oqq{x))q has the smaller cutoff effects. We expect the former to be 
the better quantity, since the latter correlator exhibits a contact term (as seen earlier), which spreads 
over a fixed number of lattice spacings. So provided the thermodynamic potentials are accurately 
known, the second strategy is the superior one. 

It remains to be seen how large the cutoff effects on J d^x{9{0) 9{x))o arc. This is shown on Fig. 5. 
For Nj- > 8, they are comparable to the cutoff effects on the T12 two-point function. Furthermore, 
they are not much larger than the cutoff effects on the entropy computed with the standard Wilson 
plaquette action [10], also displayed on the figure. 

In summary, the trace-anomaly two-point function is computationally advantageous in that one is 
computing directly a quantity which is already 0(q;^). 



3.3.2 Temperatures close to Tc 

Although perturbativc methods fail near T^, we know that very close to Tc, the specific heat becomes 
large (both in SU(3) gauge theory and full QCD), and 6s — is negative: it cancels a large flat 
contribution in the trace-anomaly correlator [11]. In that regime, it is therefore preferable to adopt 
the first strategy and compute ^ - ^ Ca^kk directly. 

3.4 Anisotropic lattice 

We now turn to the case of an anisotropic lattice, ^ > 1. Indeed it has long been recognized that such a 
lattice presents certain advantages for the calculation of thermodynamics [12] and thermal correlation 
functions, in particular in charmonium calculations [13]. 

In order to find the optimal range of anisotropies, we consider the cutoff effects on the tensor 
correlators at a fixed value of the spatial lattice spacing, a^/Lo =fixed. We then vary the anisotropy 
between 1 and 4. On Fig. 3, we see that the sign of the cutoff effects changes for 1 < ^ < 2, and goes 
to a finite value in the Hamiltonian limit, ^ ^ oo. It is clearly seen that any choice > 2 reduces 
the cutoff effects significantly. It also appears that choosing ^ > 3 does not reduce the cutoff effects 
further, presumably because they are dominated by the coarseness of the spatial discretization. The 
cutoff effects appear to be minimal near ^ = 2, and we make the choice to investigate ^ = 2 in the 
following. In fact, as Fig. (3) shows, A^^ = 16 is about as good on the ^ = 2 lattice as on the isotropic 
lattice. Because of the sign change of the cutoff effect, the smallness of the cutoff effect may be partly 
accidental. However, we find that other correlators are also improved. For instance, the discretization 
errors on the energy correlator, which are large on the isotropic lattice (Fig. 1 bottom panel), are 
significantly (Fig. 3) reduced at ^ = 2. For Nr = 20 on the anisotropic lattice, the discretization error 
is below 10% for Txq > ^, which is not the case on the Nj- = 10 isotropic lattice. We remark that 
even in the large ^ limit, the energy correlator with the present discretization is not flat, in spite of 
continuous time-translation invariance being restored in that limit. 

3.5 Non-zero spatial momentum 

For low momenta and frequencies, hydrodynamics predicts the functional form of the spectral functions 
in the shear channel (/9i3,i3) and the sound channel (p33,33) (see e.g. [14]). It is therefore of interest to 
study also correlators with non- vanishing spatial momentum [15]. An example is shown on Fig. 4 for 
p = (0,0, vrT). Here too the cutoff effects are smaller for the momentum density correlator than for 
the energy density correlator. For instance, at N^- = 16, the cutoff effects are less than 5% for Txq > ^ 
in the former case, while this level of accuracy only occurs for Txq > | in the latter case. 

4. Treelevel improvement 

Here we describe how the results obtained in this paper can be used to remove the dominant part of 
the cutoff effects on the correlators. In the case of the clover discretization, we simply divide by the 
treelevel lattice result, and multiply by the continuum result: 

C^xo, P) - a,,{xo, p) • ^ff/"'°'P^ (4.1) 

This technique is not new, see for instance [16]. 

For the plaquette discretization, the three electric-electric, magnetic-magnetic and electric-magnetic 
contributions to C(xo) are computed separately. For each of tem, we apply the following technique [9] 



to remove the tree-level discretization errors, which is adapted from static-potential studies [17]. 
Firstly, xo is defined by the treelevel correlators such that C^onti^o) = C'lat (^o)- '^^^ improved non- 
perturbative correlator C is defined at a discrete set of points through C{xo) = C(xo), and then 
augmented to a continuous function C = a + jC^^^^ . The parameters a and 7 are fixed by the condi- 
tion C(xq^) = Q! -I- 7C*ont (^0 ^ ) ) ^ = 1)2, where x^^ and x^^ correspond to two adjacent measurements. 
The correlator is thus defined for all Euclidean times between x'^^ and ^^q^^ The linear combination 
of the electric-electric, magnetic-magnetic and electric-magnetic contributions is finally obtained at a 
common value of Euclidan time. Derivatives can then also be obtained from 7^7rC'*onj(xo). 

It is clear from this discussion that it is simpler to have a site-centered definition of the energy- 
momentum tensor, and this is the choice made for our large-scale calculation [4]. 

5. Non-perturbative study of cutoff effects 

In this section we test how effective treelevel improvement is, by applying it to non-perturbative data 
from Monte-Carlo simulations. 

5.1 Isotropic lattice, plaquette discretization 

We start with data obtained on isotropic lattices with the plaquette discretization. In the scalar 
channel, the Cgg correlator was obtained at several lattice spacings in [9]. Figure 5 displays the lattice 
spacing dependence of the Txo = \ treelevel-improved tensor correlator as a function of (aT)^ for two 
fixed temperatures. The aspect ratio of the lattice is LT = 5 at 1.65rc and LT = 2.5 for 1.2ATc. The 

residual discretization errors appear to be small, in fact consistent with zero for N-j- > 8. However, a 
residual 10% cutoff effect cannot be excluded at l.GSTc, because the statistical errors are increasing 
with Nr. 

5.2 Anisotropic lattice = 2), clover discretization 

Figure 5 displays the treelevel-improved correlators of the total energy and momentum on a 16 x 48^ 
lattice. As discussed in section 3.1, these correlators are XQ-independent in the continuum limit. In 
the case of the momentum Tq^, the correlator is flat within statistical errors down to Txq = | or ^. 
Since this corresponds to xo/ua = 1.5 and 2 respectively, we regard this as an excellent outcome. By 
contrast, the unimproved correlator starts to rise for xo/ar larger by one unit. 

In the case of the energy, the treelevel improvement imposes a much larger correction to the data, 
and, not surprisingly, the expected flatness of the correlator is much less well realized. The treelevel 
improvement undercorrects the correlator for Txq > |, and overcorrects it at short distances. We have 
to conclude that in this channel, only the largest xq/ot points arc usable for a continuum extrapolation, 
even after treelevel improvement. In particular, this strongly restricts the XQ-information for finite- 
p correlators of the energy density. The latter are particularly interesting [5], because they contain 
information on the damping of sound waves in the plasma. 

6. Concluding remarks 

We have studied the discretization errors of the energy-momentum tensor correlators in lattice gauge 
theory. We summarize the lessons learnt. 

1. using a site-centered discretization of the EMT simplifies the calculation of correlators and their 
treelevel improvement. Therefore in the following we discuss the 'clover' discretization. 



2. the momentum density correlator has smah cutoff effects at treelevel, and correspondingly the 
treelevel-improved non-perturbative correlator has small cutoff effects down to very small sepa- 
rations. 

3. the energy density correlator is much more problematic with the chosen discretization. The 
treelevel energy correlator is far from being flat for realistic values of A'^,-, and correspondingly 
after treelevel improvement cutoff effects as large as 50% remain. This implies that only the 
largest values of xq can be used in a continuum extrapolation. 

4. the correlators of the spatial components of T^i, have moderate discretization errors and treelevel 
improvement works well. 

5. in all analyzed channels, the anisotropic lattice helps reduce the cutoff effects at a lower cost than 
decreasing the lattice spacing on the isotropic lattice. 

The present study suggests that the design of a lattice energy density operator which leads to 
small cutoff effects in its two-point functions would be very valuable. Such a discretization would also 
have an impact in other areas of QCD, for instance in hadron structure calculations, where the glue 
energy density operator determines the glue momentum fraction [18, 19]. It would be interesting to 
investigate the cutoff effects associated with the HYP-smeared discretizations introduced in [19]. 

A further important issue on the anisotropic lattice is the proliferation of normalization factors 
that have to be applied to different components of the energy-momentum tensor. Their determination 
will be the subject of a separate publication. 
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Figure 1: Top: the treelevel Coo.oo(a;o, 0) correlator correlator at finite lattice spacing on the isotropic lattice. 
The horizontal line is the continuum treelevel prediction, 47r^/15. Bottom: the treelevel Co3,o3(a;o, 0) correlator 
at finite lattice spacing on the isotropic lattice. The horizontal line is the continuum trelevel prediction, 47r^/45. 
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Figure 2: Top: the relative deviation of the p = tensor correlator at finite lattice spacing from the continuum 
correlator. The D's are for the clover discretization of T12 and the (g) for the clover discretization of |(Tii — T22). 
The V is the plaquette discretization. Bottom: the cutoff effects on the a^J2x(^i^)^i^))o correlator with the 
clover and plaquette discretizations. Note that C eei^N r) / Cee.cont — Coo(^T)/C'oo,cont at xq = Lq/2. The entropy 
s{Nt) is computed with the standard Wilson action and plaquette discretization of ^oo [10]. 
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Figure 3: Top: The cutoff effects on tlie tensor correlator for xq = Lq/A and Lo/2, corresponding to small 
and large symbols respectively. The spatial lattice spacing Ua is held fixed, the temporal lattice spacing is 
varied between a quarter and one times Oo-. The D's refer to T12 and the &s to ^{Tu — T22). In addition, the 
two crosses at ^ = 1 indicate the reduction of the cutoff effect on the ^(Tn — T22) correlators when increasing 
Lo/acT from 8 to 16. Bottom: the treelevel Coo,oo(a;o, 0) correlator on the anisotropic lattice, to be compared 
with Fig. (1). 
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Figure 4: The ratio of lattice to continuum treelevel correlators for Too (top) and Toa (bottom). Clover 
discretization on the anisotropic lattice, for p — {0,0, ttT). 
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Figure 5: Top: The cutoff effects on the — T22) treelevel-improved two-point function from Monte-Carlo 

simulations, for two temperatures, 1.65 and 1.24Tc. Here the plaquette discretization is used on the isotropic 
lattice. Bottom: The treelevel-improved correlators Coo,oo and J^k C'ofc.ofc at p = calculated with the clover 
discretization on the anisotropic lattice at l.STTc. These correlators are flat and equal to c„ and 3s resp. in the 
continuum limit. The unimproved correlators are displayed as open symbols. 



